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First Measurement of the ®®Ru(p, 7)®^Rh Cross Section for the p-Process with a 

Storage Ring 


Bo Thomas Aumann,^ Shawn Bishop,^ Klaus Blaum,® Konstanze Boretzky,^ Fritz Bosch,^ Carsten 

Brandau,^ Harald Brauning,^ Thomas Davinson,® Iris Dillmann,^ Christina Dimopoulou,^ Olga Ershova,^ 
Zsolt Fiilop,^ Hans Geissel,^ Jan Glorias,^ Gyorgy Gyiirky,^ Michael Heil,^ Franz Kappeler,® Aleksandra 
Kelic-Heil,^ Christopher Kozhuharov,^ Christoph Langer,® Tudi Le Bleis,^ Yuri Litvinov,^ Gavin Lotay,® 
Justyna Marganiec,^ Gottfried Miinzenberg,^ Fritz Nolden,^ Nikolaos Petridis,^ Ralf Plag,®’^ Ulrich Popp,^ 
Ganna Rastrepina,® Rene Reifarth,^’El Bjorn Riese,^ Catherine Rigollet,^® Christoph Scheidenberger,^ Haik 
Simon,^ Kerstin Sonnabend,^ Markus Steck,^ Thomas Stohlker,^’Tamas Sziics,^ Klaus Siimmerer,^ Gunter 
Weber,Helmut Weick,^ Danyal Winters,^ Natalya Winters,^ Philip Woods,® and Qiping Zhong^ 

^ GSI-Helmholtzzentrum fur Schwerionenforschung GmbH, Darmstadt, Germany 
^ Goethe University Frankfurt, Germany 
® Technische Universitdt Darmstadt, Germany 
Technische Universitdt Milnchen, Germany 
^ Max-Planck-Institut fiir Kemphysik, Heidelberg, Germany 
® University of Edinburgh, United Kingdom 
ATOMKI, Debrecen, Hungary 
^Karlsruher Institut fiir Technologic, Germany 

^National Superconducting Gyclotron Laboratory, Michigan State University, East Lansing, Michigan, USA 
KVI-GART, University of Groningen, The Netherlands 
Helmholtz-Institut Jena, Jena, Germany 

This work presents a direct measurement of the ®®Ru(p, 7)®^Rh cross section via a novel technique 
usiug a storage ring, which opens opportunities for reaction measurements on unstable nuclei. A 
proof-of-principle experiment was performed at the storage ring ESR at GSI in Darmstadt, where 
circulating ®®Ru ions interacted repeatedly with a hydrogen target. The ®®Ru(p, 7)®^Rh cross section 
between 9 and 11 MeV has been determined using two independent normalization methods. As key 
ingredients in Hauser-Feshbach calculations, the 7-ray strength function as well as the level density 
model can be pinned down with the measured (p, 7) cross section. Furthermore, the proton optical 
potential can be optimized after the uncertainties from the 7-ray strength function and the level 
density have been removed. As a result, a constrained ®®Ru(p, 7)®^Rh reaction rate over a wide 
temperature range is recommended for p-process network calculations. 

PACS numbers: 25.40.Lw, 29.20.db, 26.30.-k, 24.60.Dr 


I. INTRODUCTION 

There are 35 p-nuclei on the neutron deficient side of 
the valley of stability between ’^"‘Se and ^®®Hg, which 
are shielded against production by the neutron-capture 
processes and are produced in the p-process. The main 
production of these p-nuclei occurs via ( 7 ,n), ( 7 ,p) and 
( 7 , a) reactions, and subsequent beta decays in the so- 
called 7 -process. Therefore, it is essential to determine 
cross sections of photodisintegration reactions or their 
inverse reactions for p-process network calculations Qi- 
Thousands of nuclear reactions are involved in network 
calculations for p-process nucleosynthesis. However, 
only very few of the required cross sections have been 
measured by experiments, and thus most of them rely 
solely on predictions of the Hauser-Feshbach (HF) model 
codes, e.g., NON-SMOKER Q and TALYS 0, which 
often have very large uncertainties from nuclear input 
parameters 0 . 
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Following current p-process network calculations 
the pattern of the solar abundance for about 60% of the 
p-nuclei is reproduced within a factor of 3. However, for 
the ®®’®^Mo and ®®’®®Ru isotopes, an underproduction of 
a factor of 20-50 has been calculated for the 7 -process in 
Core Collapse Supernova models with massive stars of 13- 
25 Mq [I[. This deficiency has motivated the search for 
additional production mechanisms, e.g., the rp-process 
0 or the zzp-process Q, but also intensified efforts 
to remove the uncertainty in required nuclear physics 
parameters by measuring reaction cross sections. These 
cross sections are most sensitive to nuclear parameters, 
i.e., the 7 -ray strength function, nuclear level density 
and optical potential, in the HF model 0. However, 
these critical parameters are not well constrained by ex¬ 
periments, and hence there are large differences between 
predictions using different parameters. 

Most of the existing experimental data for the p- 
process were measured in direct kinematics using stable 
isotope targets 0. However, a direct measurement on 
unstable nuclei is still a major challenge 0 . In this work, 
we present a novel method using a heavy-ion storage ring 
developed to measure cross sections of low-energy nuclear 
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reactions, e.g., (p, 7 ) reactions, in inverse kinematics for 
nuclear astrophysics. This method offers some key ad¬ 
vantages over traditional methods. It can efficiently use 
the beam and is well suited for measurements on unstable 
nuclei with half-lives longer than several minutes. This 
method was suggested already several decades ago [1, Q . 
However, it was not realized until the latest achievements 
in producing, cooling, decelerating, and storing of heavy 
ions, as well as developments in the nuclear detection 
system at the experimental storage ring (ESR) of GSI 
[lOl| . This novel technique provides a unique condition 
for the direct measurement of (p, 7 ) reactions around the 
energy range of astrophysical interest m and has been 
successfully demonstrated for the first time by measuring 
the ®®Ru(p, 7 )®’^Rh cross section between 9 MeV/u and 
11 MeV/u. A preliminary data analysis at 11 MeV/u 
has been reported in conference proceedings [ 13 , fl3| . 
Here we report a full analysis with all corrections made. 
Experimental results of the present work allow us to 
constrain the most important parameters in the HE 
model, and thus provide a reliable prediction for this 
reaction over a wide energy range. 



FIG. 1 . (Color online) The experimental setup for proton 
capture reactions in inverse kinematics at the ESR. The H2 
target and the detectors used in this experiment are marked. 
Both the Multi-Wire Proportional Chamber (MWPC) and 
the double-sided silicon strip detectors (DSSSDs) were placed 
in pockets with 25 /rm thick windows. A schematic view of the 
DSSSDs is also shown. The unreacted ions (magenta 

solid line after the target) were recycled and focused on the 
H2 target repeatedly during the measurement phase. 


II. EXPERIMENTAL DETAILS 

During this experiment, ®®Ru ions from the linear ac¬ 
celerator (UNILAC) were first accelerated to 100 MeV/u 
in the heavy-ion synchrotron (SIS) and then stripped to 
the bare charge state of 44-1- using all mg/cm^ carbon 
stripper foil. The fully stripped ions were injected into 
the ESR and slowed down to 9 MeV/u, 10 MeV/u and 
11 MeV/u, respectively, by ramping the magnetic fields 
and the frequency of the radio-frequency (rf) system 
synchronously. However, large beam losses occurred 
during this deceleration phase mainly due to imperfec¬ 
tions of the ramping parameters. The ®®Ru'^^“'" ions 
were cooled by the electron cooler before and after the 
slowing down phase to a small diameter (about 5 mm) 
and momentum spread (around 10“^). After the final 
slowing down phase, about 5x10® ®®Ru‘^^“'' ions were 
stored in the ESR with a lifetime of several hundred 
seconds when the hydrogen target was switched off. 
Finally, a windowless hydrogen microdroplet target of 
high density [l 3 | was switched on and the decelerated 

96Ru44-e 

ions were focused onto this target for nuclear 
reactions. The great advantage of this storage ring 
method is that unreacted ®®Ru^"'^“'' ions were recycled and 
repeatedly impinged on the hydrogen target for reactions. 
Considering the revolution frequency of about 400 kHz 
for ®®Ru"'^^“'' ions circulating in the ring and the thickness 
of the H 2 target of about 10 ^® particles/cm^, a luminosity 
of about 2 x 10 ^® /cm^/s has been achieved. 

At about 10 MeV/u, the main reaction channels of 
^^Ru^^"*" with the hydrogen target include the atomic 
electron capture (EC) reactions and different nuclear 
reactions. The former contain mainly two parts, namely 
the non-radiative electron capture (NRC) and the radia¬ 


tive electron capture (REC) accompanied by the emission 
of a photon. There are large uncertainties of about 
30% according to Ref. [l^ in accurately determining 
the absolute beam intensities, target densities and the 
beam-target overlap at the ESR. To remove these large 
uncertainties, both the K-shell REC (K-REC) products 
(photons) and the EC products were registered by our 
detectors, which allowed us to absolutely determine the 
(p, 7 ) cross section by using two normalization methods. 

Fig.[T]shows the experimental setup from the hydrogen 
target to all used detectors at the ESR. The X-rays 
emitted from the atomic reactions were registered by 
a Ce detector mounted close to the target interaction 
area at the observation angle of 90° with respect to the 
beam axis. The detection efficiency of this Ce detector 
was calibrated with mixed 7 -ray sources in the energy 
range between 10 keV and 140 keV. In the K-REC 
energy region, an intrinsic efficiency of about 88 % has 
been reached for this detector with a solid angle (AH) 
of about 2.5x10“® sr. The down-charged ®®Ru ions 
(96Ru43-e, 96Ru42-h^ 96Ru4i+, etc.) produced by EC were 
recorded by a position sensitive Multi-Wire Proportional 
Chamber (MWPC) [l^l mounted in the vacuum chamber 
of the first dipole magnet behind the target, as indicated 
in Fig. [H This detector was operated with a mixture 
of argon, CO 2 and heptane (80:20:1.5) gas at standard 
atmospheric pressure in a pocket with a stainless-steel 
window. The detection efficiency of this detector is better 
than 99% for ions above ^10 MeV/u and a position 
resolution (FWHM) of 1.9 mm has been reached M- 
However, a fraction of electron capture events were lost 
at the beam energy of 9 MeV/u due to significant energy 
losses in the pocket window with a thickness of about 
25 pm and the gas with a thickness of around 24 mm 
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before the MWPC. 

At about 10 MeV/u, there are only four open nu¬ 
clear reaction channels: ®®Ru(p,p)®®Ru elastic scatter¬ 
ing, 96 Ru(p, 7 ) 97 Rh, 96 Ru(p, a)93Tc and ^^Ru{p, nf^Rh. 
The above reaction products with different mass-over¬ 
charge ratios (m/q) were separated by the magnetic field 
{B) of the dipole magnet behind the target. For instance, 
the dispersion of the magnets displaced ®^Rh^3+ from 
unreacted ®®Ru'^'*“'- by about 142 mm at the position of 
our detector, see dash-dotted line and solid line in Fig.[T] 
As will be shown below, the (p, 7 ) reaction products can 
be discriminated from other nuclear reaction products 
due to their relatively small momentum spread. The m/q 
is always larger for atomic EC reaction products since 
nuclear reaction products are bare. Therefore, orbits of 
the former are always on the outer side of the ESR, see 
dashed line in Eig. [TJ and bare nuclear reaction products 
on the inner side are not contaminated with atomic ones. 

The separated nuclear reaction products were detected 
by the position sensitive double-sided silicon strip detec¬ 
tors (DSSSDs) behind the quadrupole triplet magnets, 
as shown in Fig. [TJ The DSSSDs consist of two silicon 
detectors with a ^^ 21 .5 mm inactive gap between them, 
see Fig. [TJ Each silicon detector with an active area of 
50x50 mm^ has 16 strips in both x- and j/-directions. 
The strip pitch is 3.1 mm and the strip length is 49.5 mm. 
The DSSSDs were placed in a pocket separated from 
the vacuum of the ESR by a 122x44 mm^ stainless- 
steel window with a thickness of about 25 /rm. The 
thickness of this window limited the reaction energy to 
above 9 MeV/u in this experiment, since heavy ions with 
less energy would already be stopped in the window. 
Recently, an improved detector omitting the window has 
been mounted [l7l| , which will allow us to measure nuclear 
reactions around 5 MeV/u in future experiments. 


III. DATA ANALYSIS AND RESULTS 
A. Experimental spectra 

Eor each beam energy setting, the DSSSDs were moved 
to two different positions along the x direction, e.g., 
0 mm and 25 mm, to combine two measured spectra 
into a common x position spectrum without any gap, 
see black points at 11 MeV/u in Eig. |2J To identify 
the (p, 7 ) reaction products unambiguously and to study 
their transmission efficiency, a Geant4 [l8l| simulation 
has been performed using a numerical model of the 
experimental setup shown in Eig. |T| Fig. [5] compares the 
simulated total x position spectrum including all nuclear 
reaction events (red line) with the experimental data 
(black points) registered by the DSSSDs at 11 MeV/u. 
The experimental data, which have an uncertainty of 
about 10 %, can be well reproduced and the (p, 7 ) reaction 
events can be disentangled clearly from the background 
produced by (p,p), (p, o), and (p, n) reactions, based on 
Geant4 simulations. The angular distribution of (p,p) 
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FIG. 2. (Color online) The experimental x position distribu¬ 
tion of nuclear reaction products (black points) registered by 
the DSSSDs at 11 MeV/u is compared with the simulated x 
position distribution (red line). Events from different reaction 
channels, i.e., (p,p), (p, a), (p, 7 ), and (p, n) reactions, can 
be disentangled based on the Geant4 simulation. The (p, 7 ) 
reaction products have a narrower distribution than other 
nuclear reaction products due to their smaller momentum 
spread. 


Experiment at 11 MeV/u 

- Simulation at 11 MeV/u 

- Simulated for (p,p) 

- Simulated for (p,n) 

- Simulated for (p,aj 
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scattering products in the center-of-mass (CM) system, 
which serves as an input for the simulation, is from 
the prediction by TALYS [l^ while other events are 
assumed to be isotropic. The magnetic fields in the 
simulation were set to the experimental values. The 
minimum chi-square method has been utilized to obtain 
the best simulation spectrum for the experimental data 
and determine the number of (p, 7 ) products The 

simulation uncertainty can be obtained by simulations 
varying the sensitive parameters, e.g., magnetic fields, 
sizes of beam pipes and sizes of chambers, within their 
uncertainties. 


The charge-state spectra measur ed by the MWPC at 
11 MeV/u are presented in Figure 3(a) Five different 
charge states of ®®Ru are produced by the electron 
capture reactions. The most prominent peak is from 
93 Ru 43+ ions produced by single EC. Other peaks are 
caused by the capture of more electrons. The black 
line shows the spectrum measured when the H 2 target is 
switched on while the red line indicates the background 
measured when the H 2 target is off. Hence, the single 
EC events caused by the H 2 target can be determined 
by subtracting the corresponding background from the 
interaction of the ^^Ru^^"*- beam with the residual gas in 
the ring. 


The X-ray spectrum measured by th e 90° Ge detector 
at 11 MeV/u is given in Figure [3(b)[ The K-, L- and 
M-REC peaks are caused by the radiative captures into 
the K-, L- and M-shells of Ru, respectively. The Kq,, 
Rp and K.^ peaks originate from cascades after electron 
captures into higher shells of Ru. Positions of these peaks 
are in very good agreement with theoretical predictions 
of Ref. (20|. According to our experimental data with 
the H 2 target switched off, the background spectrum 
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FIG. 3. (Color online) (a) Position spectra of EC events 
measured by the MWPC at 11 MeV/u. The black and red 
lines refer to the spectrum with and without the hydrogen 
target, respectively. Spectra have been downscaled by a factor 
of about 300. (b) X-ray spectrum (black line) at 11 MeV/u 
registered by the 90° Ge detector. The spectrum is not 
corrected for the detection efficiency. The background is fitted 
with a linear function. 


registered by the 90° Ge detector is almost linear and 
there is no peak structure in the K-REC energy region. 
Therefore, the number of K-REC events induced by 
the H 2 target can be extracted from the sum of all 
events in the K-REC energy region, subtracting a linear 
background. 


B. Cross section determination 


In this experiment, the (p, 7 ) cross section can be 
normalized by two methods using ( 1 ) the theoretical 
single EC cross section, and (2) the calculated K-REC 
cross section at 90°. The single EC cross section and 
the K-REC cross section can be predicted very well by 
different theoretical models. 

For the single EC at about 10 MeV/u, the cross section 
can be calculated by the Schlachter scaling rule [2l|. 
According to this scaling rule, the cross section of the 
single EC can be calculated by the relation 


CTEC = 1-1 X 


-4.8 


barn. 


( 1 ) 


where q is the projectile charge state, Z 2 is 1 for the 
hydrogen target, and E is the projectile energy in keV/u. 


TABLE I. Theoretical cross sections of the single EC and 
corrected differential K-REC cross sections at 90° for 
-I- H 2 collisions between 9 MeV/u and 11 MeV/u. Single EC 
cross sections are also normalized to corrected K-REC cross 
sections at 10 MeV/u and 11 MeV/u. 


E (MeV/u) CTEC (barn) duK-REc/dll (barn/sr) (Teq™ (barn) 


9 

2960±592 

67.06±13.42 


10 

1780±356 

59.72±11.95 

1403±351 

11 

1130±226 

53.74±10.75 

1272±318 


The calculated single EC cross section for our experiment 
is given in Table ID Its uncertainty is estimated to be 
about 20%. The single EC cross section above 10 MeV/u 
can also be normalized to the theoretical K-REC cross 
section, see below, by a validated normalization method 
reported in Ref. [1^ when both single EC and K-REC 
events were recorded by our detectors. Using this nor¬ 
malization method, the single EC cross section has been 
calculated by Eq. (2.2) and Eq. (2.3) in Ref. [l^ (see also 
Eq. ([5]) given below with the (p, 7 ) replaced by the single 
EC). Single EC cross sections determined by two different 
methods are in good agreement within uncertainties, see 
Table H This agreement also indicates the Schlachter 
scaling rule works well for ^^Ru"^^"*" colliding with H 2 at 
about 10 MeV/u in our experiment. 

With the number of the (p, 7 ) reaction events 
and the number of the single EC events Vec, the (p, 7 ) 
cross section can be normalized to the single EC cross 
section ctec, which has been calculated by the scaling 
rule, using 




( 2 ) 


and Nec can be derived from Fig.[5]and Fig. 3(a)[ 
respectively. However, the (p, 7 ) cross section at 
9 MeV/u cannot be determined by this method because 
some electron capture events produced at this energy 
were stopped in the pocket window and gas in the pocket 
before reaching the detector. 

Since the K-REC cross section is well-known, the 
number of K-REC events registered by the 90° Ge 
detector is also used to accurately calculate the (p, 7 ) 
cross section by means of the normalization method. The 
REC is the inverse of the photoelectric effect, which can 
be described in the framework of a nonrelativistic dipole 
approximation by Stobbe [l^. Thus, the total K-REC 
cross section can be obtained by the principle of detailed 
balance according to [l^ 


ck-rec = 9165 X 


\ exp[—4i/arctan(l/j/)] 
1 + J 1 — exp( —2711/) 


(3) 

1 / = aZjj3 is the Sommerfeld parameter, a = 1/137.036 
is the fine-structure constant, f3 = v/c is the projectile 
velocity, and Z is the projectile atomic number. Within 
the dipole approximation, the differential K-REC cross 
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section at 90° can be calculated by the formula [T^ 

dcTK-REC^a n^,o^ ^ 

_^(0 = 9O ) = -aK-REC, (4) 

where 0 is the X-ray observation angle with respect to the 
beam direction. Then one can calculate the differential 
K-REC cross section at 90° by Eq. ([3]) and Eq. (jH). 
However, it has been found that the theoretical value 
predicted by this formula overestimates the experimental 
K-REC cross section by a factor of about 30% MM, 
which appears to be a common feature for K-REC 1^ . 
Therefore, an empirical factor of 0.7 for the theoretical 
prediction was recommended by Stohlker et al. [H, [2^ . 
Theoretical cross sections adjusted by this factor agree 
very well with measured K-REC cross sections, see, e.g., 
Ref. [13. Thus, this factor has also been applied in our 
calculations of K-REC cross sections. The differential K- 
REC cross section at 90° corrected by this factor is listed 
in Table m 

With iV(p,-y) and the number of K-REC events regis- 
tere d by the 90° Ge detector A^k-rec, extracted from 
Fig. |3(b)| the (p, 7 ) cross section can be normalized to 
the differential K-REC cross section at 90° duK-REC /dH 
by the following expression 



Energy (MeV) 

FIG. 4. (Color online) A comparison of the measured {p, 7 ) 
cross sections with the predictions by the standard NON- 
SMOKER code [3 (dash-dotted line) as well as the TALYS- 
1.4 code [3 using different 7 -ray strength functions, i.e., the 
HFB (solid line) [13, Brink-Axel Lorentzian (dotted line) 

a [13 and Kopecky-Uhl generalized Lorentzian (dashed line) 
. Red points and purple squares represent the measured 
{p, 7 ) cross sections normalized to the K-REC cross sections 
and the single EC cross sections, respectively. The latter 
(purple squares) have been slightly offset in energy for clarity. 
The inset shows a zoom around our measured data. Note that 
all energies are in the CM system. 


^(p.7) 




^(P.7) „ dcTK-REC,^^ 

—EK-REC- 


K-REC 


dil 


(5) 


Ek-rec is file intrinsic efficiency and AO is the solid angle 
spanned by the Ge detector. The K-REC cross section 
listed in Table |T] is adopted to calculate the {p,j) cross 
section using Eq. ©■ Similar methods have been applied 
in atomic K-REC experiments at the ESR and proved to 
be justified [T^ - 


The (p, 7 ) cross sections normalized to K-REC cross 
sections (red points) are in good agreement with those 
normalized to single EC cross sections (purple squares), 
as displayed in Fig. ID This agreement gives confidence in 
the present normalization methods. The uncertainty of 
the (p, 7 ) cross section is dominated by the uncertainty 
of the calculated K-REC or EC cross section (20%) and 
the uncertainty of -j,) (^ 20 %) while the uncertainty 
of other factors is always less than 10%. In the following, 
(p, 7 ) cross sections normalized to K-REC cross sections 
will be compared with HE calculations since K-REC cross 
sections have been checked by many atomic experiments 
and proved to agree well with the experimental values 
within a small uncertainty of about 20% [l3[23- Another 
reason is that K-REC events have less background than 
EC events, see Figs. 3(a) and |3(b)[ The ®®Ru(p, 7 )®^Rh 
cross section normalized to the K-REC and the associ¬ 
ated S factor are summarized in Table [Hi 


TABLE 11. Measured cross sections and astrophysical S 
factors for ®®Ru(p, 7 )'^^Rh. 

Ecm (MeV) Cross section (mbarn) S factor (10® keV barn) 


8.976 

Q 00+2.58 

0.Z5_2.76 

i.49t°:f 

9.973 

7 00+2.13 

' •oO_2.i3 

0.74t°-2 

10.971 

Q 10+2.59 
^•^'^-2.94 



IV. MODEL CALCULATIONS 

A. Cross sections 

The ®®Ru(p, 7 )®^Rh cross section is mainly sensitive to 
the 7 -ray strength function in our experimental energy 
region (^10 MeV in the CM system) while it depends 
sensitively on both the 7 -ray strength function and the 
proton potential in the Gamow window between rougUw 
1.3 MeV and 4.3 MeV, see e.g.. Fig. 16 in Ref. 0. 
Therefore, our experimental results can be used directly 
to remove the uncertainty from the 7 -ray strength func¬ 
tion and improve the agreement between the theoretical 
predictions and the experimental data over a large energy 
region. The comparison of our experimental results by 
two normalization methods with predictions by NON- 
SMOKER [3 using default parameters as well as predic¬ 
tions by the TALYS-1.4 code [3 utilizing different 7 -ray 
strength functions is presented in Fig. |4l 

The default proton potential from the parameteriza- 
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tions of Koning and Delaroche (KD) [2^ is adopted 
in TALYS using different El 7 -ray strength functions, 

1. e., the Kopecky-Uhl generalized Lorentzian (default 

strength) [ 23 , Brink-Axel Lorentzian (BAL) and 

microscopic Hartree-Fock-Bogoliubov (HFB) [^. As 
shown in Fig.lH the (p, 7 ) cross section is sensitive to the 
7 -ray strength function from about 2 MeV to very high 
energies and differences between cross sections predicted 
by TALYS using different 7 -ray strength functions are 
particularly large at energies above 5 MeV. Therefore, 
our experimental results between 9 MeV and 11 MeV 
provide an excellent constraint for the 7 -ray strength 
function, which can also have a significant impact on the 
(p, 7 ) cross sections in the Gamow window. 

The microscopic optical potential of Jeukenne et al. 
[ 2 ^ is applied in the standard NON-SMOKER code. The 
7 -ray strength function and nuclear level density models 
in NON-SMOKER are from Ref. [s^ and Ref. [3l|, 
respectively. The (p, 7 ) cross section predicted by NON- 
SMOKER using the above input has a typical uncertainty 
of a factor of 2 [s^. Predictions by this code agree well 
with our experimental results at 9 MeV and 10 MeV. 
However, NON-SMOKER shows a very drastic tendency 
of decreasing from 9.5 MeV to higher energies and thus 
underestimates our data at 11 MeV by a factor of about 

2. The possible reason for this discrepancy is that the 
7 -ray width is artificially downscaled at high energies to 
simulate the effect of the pre-equilibrium and in-cascade 
particle emission in this code [33|. 

For our experimental results, the best agreement is 
achieved by TALYS using the BAL 7 -ray strength func¬ 
tion. Other models implemented in TALYS severely 
underestimate our results. A slight deviation between 
predictions by TALYS using the BAL strength function 
and our experimental results at 11 MeV may stem from 
the uncertainty of the proton potential, see below. Con¬ 
sidering the good agreement, the BAL strength function 
will also be applied to predict the (p, 7 ) cross sections 
at lower energies. Note that the uncertainty from the 
neutron potential has been checked by varying the default 
KD neutron potential by a factor of ^100 since the (p, 7 ) 
cross section shows a small sensitivity to it. It is found 
that the variation of the (p, 7 ) cross section is very small 
and thus the uncertainty from the neutron potential does 
not change the above conclusions. In Fig. m the nuclear 
level density adopted in TALYS calculations is the widely 
used back-shifted Fermi gas (BSFG) model (s^, [s^. 

The (p, 7 ) cross section also shows a relatively small 
sensitivity to the nuclear level density model. Thus, 
different level density models, i.e., the Constant Tem¬ 
perature (CT) [ 3 ^, BSFG [SJ, 113 , Generalized Super¬ 
fluid (GS) [131 and Hartree-Fock-Bogolyubov (HFB) plus 
combinatorial [s^ models, have been tested in TALYS 
by using the same BAL strength function, as shown in 
Fig. [ 5 ] Both the BSFG model and the GS model agree 
rather well with our experimental data. The relative 
difference between these two models is always less than 
5%. Therefore, Fig. [S] also supports our conclusion that 



FIG. 5. (Color online) A comparison of the measured 
(p, 7 ) cross sections with the predictions by the TALYS- 
1.4 code [3 using different level density models, i.e., the 
Constant Temperature (CT) [ 3 ^, BSFG [ll. [35|. Generalized 
Superfluid (GS) [37|and Hartree-Fock-Bogolyubov (HFB) 
plus combinatorial models. The BAL 7 -ray strength 
function is applied in the TALYS calculations. The purple 
squares have been slightly offset in energy. 

TALYS using the BAL 7 -ray strength function as well 
as the BSFG (or GS) level density model gives the best 
agreement with our experimental data between 9 MeV 
and 11 MeV. 

After the 7 -ray strength function as well as the level 
density have been strongly constrained by our experimen¬ 
tal data for ®®Ru(p, 7 )®’^Rh in the previously unexplored 
energy region, some additional data for this reaction 
between 1.6 MeV and 3.4 MeV measured by J. Bork 
et al. via the activation method can help to improve 
the prediction of the ®®Ru(p, 7 )®^Rh cross section further 
into the Gamow window, where it is sensitive to both 
the 7 -ray strength function and the proton potential [ 3 ] . 
Recent experimental studies for other proton induced 
reactions indicate that the proton potential should be 
modified to explain the measured cross sections [40l - l43| . 
The reason is that the parameters used in the proton 
potential are usually derived from reactions at energies 
far above the astrophysical energy region [3|. In the 
following, the modification of the proton potential will be 
tested and the uncertainty from the proton potential will 
be constrained by comparing with experimental results 
for 96Ru(p,7)97Rh. 

In Fig. E 9®Ru(p, 7 ) 97 Rh cross sections from two 
experiments in different energy regions are compared 
with predictions by NON-SMOKER Q and TALYS Q 
using the standard KD proton potential as well as TALYS 
using the modified KD proton potential. In the modified 
KD proton potential, default values of ry and di have 
been scaled by factors of 0.7 and 0.8, respectively, to 
obtain a good agreement with the data measured by 
Bork et al. [13] and also our new data. The BAL 7 -ray 
strength function and the BSFG level density are utilized 
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Energy (MeV) 

FIG. 6. (Color online) Our experimental results (red points) 
and the experimental results of J. Bork et al. (black points) 
[ 3 ^ in different ener^ regions are compared with predictions 
by NON-SMOKER (blue dashed line) and TALYS using 
the KD proton potential (black dash-dotted line) as well as 
TALYS using the modified KD proton potential (magenta 
solid line). The staggering structure in the low energy 
experimental data is due to the low level density. Note 
that both the cross section and the energy are plotted on a 
logarithmic scale. 


in TALYS because they can provide the best prediction 
for our results between 9 MeV and 11 MeV, see Figs. |4] 
and El Both NON-SMOKER and TALYS using the 
standard KD proton potential significantly overestimate 
cross sections at low energies, e.g., by a factor of around 
5 at energies below 2 MeV. TALYS using the BAL 
7 -ray strength function, the BSFG level density and 
modified KD proton potential can excellently reproduce 
both experimental results in different energy regions. 
It is remarkable to see that TALYS using constrained 
parameters can reproduce cross sections varying from 
about 3x10“® barn to ^1x10“^ barn in a large energy 
range from about 1.5 MeV to 11 MeV. 


B. Reaction rates 

As an input for p-process network calculations, the 
stellar rate Na{<J* v )over a large temperature range, 
especially in the astrophysically relevant temperature 
range between 1.5 GK and 3.5 GK, is required [H, 
[^. However, the ®®Ru(p, 7 )®^Rh cross section between 
1.6 MeV and 3.4 MeV measured by J. Bork et al. is 
only sufficient to compute the rate between about 2 GK 
and 2.5 GK [s^, see black triangles in the filled region in 
Fig.m Experimental data must be extended to a wider 
energy region by HF calculations. The TALYS cross 
section constrained by two experimental data sets, as 
presented in the previous section, represents a solid basis 
for the reaction rate calculation covering the important 



Temperature (GK) 

FIG. 7. (Golor online) Stellar reaction rates over a large 
temperature region calculated by the TALYS code using 
the modified KD proton potential (red solid line). They 
are compared to rates from the experimental data of Bork 
et al. (black triangles in the filled region) as well as 

the theoretical extrapolation (black triangles out of the filled 
region), NON-SMOKER calculations (blue dashed line) 
and BRUSLIB (magenta dash-dotted line). Both the rate 
and the temperature are illustrated on a logarithmic scale. 


TABLE III. Recommended stellar rates at different temper¬ 
atures for ®®Ru(p, 7 )®^Rh. The rate has an uncertainty of 
about 30% according to the comparison with measured data. 


Temperature (GK) 

Stellar rate (cm^s ^mole ^) 

0.5 

(1.8±0.5)xl0““ 

0.6 

(9.1±2.7)xl0“i'’ 

0.7 

(2.0±0.6)xl0“® 

0.8 

(2.5±0.8)xl0“^ 

0.9 

(2.1±0.6)xl0“® 

1 

(1.3±0.4)xl0“® 

1.5 

(7.0±2.1)xl0“® 

2 

(3.4±l)xl0“^ 

2.5 

5.0±1.5 

3 

(3.7±l.l)xl0^ 

3.5 

(1.8±0.5)xl0^ 

4 

(6.1±1.8)xl0^ 

5 

(3.7±l.l)xl0® 

6 

(l.l±0.3)xl0'‘ 

7 

(2.3±0.7)xl0'‘ 

8 

(3.4±l)xl0'‘ 

9 

(4.1±1.2)xl0'‘ 

10 

(4.4±1.3)xl0'‘ 


temperature range. 

The stellar reaction rates calculated by the TALYS 
code using constrained nuclear physics input are com¬ 
pared to the experimental rates around 2 GK from 
Ref. as well as their theoretical extrapolation, NON- 
SMOKER calculations Q and BRUSLIB [i^, as pre¬ 
sented in Fig.jT] The experimental rates have been theo- 










TABLE IV. Recommended REACLIB parameters for the reactivity of ®®Ru(p, 7 )®^Rh and its reverse reaction. 



ao 

fli 

a2 

as 

04 

as 

as 

(P,7) 

0.90232 

-10.08917 

-5.75549 

2.62112 

1.01194 

-0.29641 

4.64569 

(7>P) 

22.29232 

-54.302315 

-5.75549 

2.62112 

1.01194 

-0.29641 

6.14569 


retically extrapolated by normalizing the ®®Ru(p, 7 )®^Rh 
cross sections calculated by NON-SMOKER to the ex¬ 
perimental data of J. Bork et al., see Ref. [s^. In 
the normalization, a factor of about 0.5 has been ap¬ 
plied for NON-SMOKER calculations. The TALYS 
rate using parameters constrained in the previous sec¬ 
tion can excellently reproduce the experimental rate for 
9®Ru(p, 7 ) 9 ^Rh between 2 GK and 2.5 GK, as displayed 
in Fig. [71 Hence, the TALYS rate constrained by two 
experimental data sets at different energy regions is 
recommended for ®®Ru(p, 7 )®^Rh, as listed in Table m 

On the contrary, both NON-SMOKER and BRUSLIB 
overestimate the rate at temperatures below 3 GK, 
particularly at low temperatures. A good agreement 
is reached between predictions by NON-SMOKER and 
TALYS between 3.5 GK and 9 GK. However, NON- 
SMOKER underestimates the rate above 9 GK, which 
is caused by the underestimation of the cross section 
above about 10 MeV. Above 3 GK, the rate extrapolated 
by NON-SMOKER is lower than the TALYS rate since 
the former are determined by normalizing the NON- 
SMOKER rate to the experimental rate around 2 GK. In 
addition, BRUSLIB significantly underestimates the rate 
above 4 GK. For instance, BRUSLIB underestimates the 
rate by a factor of 20 above 8 GK, which again indicates 
that our experimental data are important to constrain 
the theoretical rate. It should be stressed that the 
BRUSLIB rate is calculated by the TALYS code, where 
the input parameters are not constrained by experiments, 
see Ref. for details. 

For network calculations, the recommended rate has 
been parameterized in the REAGLIB format Q using 
the formula 

Na{(Jv) = exp[ao -I- aiTg“^ -|- 

+a 4 Tg + a^Tg^ -I-og In (Tg)], ( 6 ) 

where Tg is the temperature in GK. Recommended REA¬ 
CLIB parameters are listed in Table lYl The fit using 


these parameters agrees very well with the recommended 
rate within 10% between 1 GK and 8 GK. 


V. SUMMARY AND OUTLOOK 

In summary, a novel technique via the collision of 
stored heavy ions with a hydrogen target has been devel¬ 
oped at the ESR, which provides unrivalled opportunities 
for the direct measurement of (p, 7 ) reactions around the 
energy range of astrophysical interest, particularly for 
previously unreachable radioactive ions. This method 
has been successfully demonstrated for the first time 
by measuring the ®®Ru(p, 7 )®^Rh cross sections between 
9 MeV and 11 MeV. The present experimental results 
allowed us to pin down the 7 -ray strength function, 
which is a critical parameter in the HE model, as well 
as the nuclear level density model. After this, another 
important parameter, the proton potential, has also been 
constrained by combining our results with some addi¬ 
tional data at lower energies. TALYS, constrained by two 
experiments in different energy regions, can excellently 
predict the stellar rates for ®®Ru(p, 7 )®^Rh over a large 
temperature range for p-process network calculations. 

Further measurements of (p, 7 ) reactions at lower en¬ 
ergies around the Gamow window via this method using 
our improved detector are planned in future experiments 
at heavy ion storage rings. Besides, ( 01 , 7 ) reactions will 
also be measured by this method when a helium target 
is utilized. 
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